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Abstract 

We prove sharp rates of convergence to the Ewens equilibrium distribution for 
a family of Metropolis algorithm based on the random transposition shuffle on the 
symmetric group, with starting point at the identity or at a random n-cycle. The 
proofs rely heavily on the theory of symmetric Jack polynomials, developed initially 
by Jack [T5], Macdonald [19] . and Stanley [22]. This completes the analysis started by 
Diaconis and Hanlon in [5]. In the end we also explore other integrable Markov chains 
that can be obtained from symmetric function theory. 
> ■ 

: 1 Introduction 

The connection between the symmetric group S n and the set of integer partitions of n is 
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widely known. The partition that a permutation a G S n corresponds to is simply given 
by its cycle structure. In fact this connection is the basis for the classical representation 
theory of S n (see e.g. [10]): the set of irreducible representations of S n are indexed by 
the set V n of partitions of n. Since S n is finite, it can also be endowed with a probability 
space structure. The most natural measure on S n is thus the uniform measure, with each 
permutation getting a weight of l/n\. Sampling from this uniform measure is important for 
many statistical applications ([!]), such as testing independence of n iid uniform random 
variables on an ordered set. Its intimate connection with card shuffling models has also 
generated a wonderful array of mathematics, most notably the determination of their mixing 
times (see for instance [2], [8], and [21]). 

One of the most natural generalizations of the uniform measure on S n is a 1-parameter 
family of so-called multivariate Ewens distributions, named after Warren Ewens, who derived 
the partition function of this probability measure. It is defined by giving each permutation 
a a weight of a l ^\ a > 0, where £(a) is the number of cycles in a; hence it can be viewed 
as an exponentially tilted family based on the uniform measure. The distribution was first 
applied to population genetics, in which it describes the distribution of frequencies of alleles 
in a sample of genes ([TB] Chapter 41). 
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Many important properties of the uniform measure on S n continues to hold for the 
multivariate Ewens distribution ([I]). For instance, the two perfect sampling schemes of the 
uniform measure, Feller coupling and Chinese restaurant process, both admit generalizations 
to the Ewens case. In this article, we describe a much more subtle property of the uniform 
measure that has been successfully generalized to the a deformed setting. In a nutshell, the 
characters of the symmetric group S n form a basis in the Fourier space of class functions on 
S n under the uniform measure [7J. When the underlying measure is a-deformed from the 
uniform, one can only make sense of Fourier transforms of a particular type of class functions, 
namely the ones supported on transpositions and the identity class. In that case, the basis 
in the Fourier space becomes the matrix coefficients of the transition from Jack symmetric 
polynomials basis to the power sum symmetric polynomials basis. These generalize the 
characters, which happen to be the transition coefficients from Schur polynomial basis to 
power-sum basis. 

This Fourier analytic property was first obtained by Stanley in [22]. Later Hanlon (|12j) 
applied it to the study of the Metropolis Markov chain based on random transposition walk on 
Syr t licit converges to the multivariate Ewens distribution. Diaconis and Hanlon ([5]) further 
initiated the investigation of total variation mixing time of this chain. They conjectured the 
upper bound to be of the same order as in the uniform case, which we prove rigorously here. 
Also we provide a matching lower bound that together with the upper bound implies the 
cut-off phenomenon for the a-deformed chain, first conjectured in [5] as well. 

Next we study the case of starting the Metropolis chain at the other extreme, the uni- 
form distribution on the set of n-cycles. The behavior is patently different from starting at 
the identity: the mixing time is of order 0(n), and there is no cut-off in the distance to 
stationarity profile. 

In the appendix, we include some preliminary attempts to generalize the walk studied 
here in various directions. First we look at the action of the Sekiguchi-Debiard operator on 
other classical bases of symmetric polynomials. We also consider higher order operators, as 
given by the operator valued generating function in [19] page 317. These turn out to give 
new local move Markov chains converging to MED (9), albeit without simple group theoretic 
interpretations. Lastly we look at Laplace-Beltrami operators associated with other root 
systems. Recall the Schur- Weyl duality between the simple Lie groups SU (n) and the finite 
groups S n . This leads to an interpretation of the Sekiguchi-Debiard operator (as well as 
their quantized version given by Macdonald [12] ) as associated with root system of type A n . 
The appropriate generalizations were first discovered by Heckman and Opdam ([13]) in the 
context of Hamiltonian systems of particles on a circle, and later extended to the Macdonald 
case in [20] (see also [17] and references therein for a 5-parameter generalization). 

Jack polynomials, which form the backbone of the argument presented here, turn out to 
be special cases of Macdonald polynomials of type A n . Diaconis and Ram [B] interpreted 
them as eigenf unctions of an auxilliary variable algorithm on the space of partitions, which 
can be viewed as a quantized version of the local walk studied here. 
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2 Metropolis walk starting at the identity class 



The multivariate Ewens distribution with parameter a is defined on S n with P(o~) pro- 
portional to o/^ where £(cr) is the number of cycles of a. The normalization constant 
z n (a) = ac(n) '■= a{a + 1) . . . (a + n — 1). 

Consider now the random transposition walk on S n , defined by picking a pair of numbers 
i 7^ j at random, and multiplying the current state in S n by the transposition (ij). In this 
form, the walk is periodic and does not converge. But if we make it lazy, then it converges to 
the uniform measure on S n . By metropolizing the non-lazy walk to the multivariate Ewens 
distribution with parameter a, we create a new markov chain that converges to MED(a) 
(see [5] for details). The walk behaves differently for a > 1 and a < 1. And as long asa^ 1, 
it converges, because it always has positive holding probability. 

Warning. The a parameter here will be the reciprocal of 9 below. 

If we start the Metropolis walk described in the introduction from the identity element, 
then we can view it as a walk either on the symmetric group or on the set of partitions. It 
is the latter interpretation that allows for sharp analysis with other starting points. 

Theorem 2.1. For 9 G (0, oo) \ {1}, let Pg be the discrete time Metropolis chain based on 
random transposition walk starting from id converging to the multivariate Ewens distribution 
it with parameter 9" 1 (so identity has the largest mass when 9 < 1). Explicitly, let X(it) be 
the cycle structure of the permutation n, X' be the tranposition of X as a Ferrers diagram, 
and denote n(X) = Y^=i (dO ■ V ^ e t rans tt>ion rule is given by 

' l-lA^ + ^l^lAfl-lAO ifo = Tx 

7sy(lA0) 2 ifa = 7t(i,j) and 1 

|!j(lAn if£(a)=£(n) + l 

otherwise 



it, a 



then the chain has a cutoff at t = |(| V l)nlogn. This means 

lim limsup H-f^ — 7r|| TV = 



lim liminf ||Py — 7t||tv = 1 

c— >— 00 n 

for t(c) := |(| V l)n(logn + c). 

Remark 1. 1. Notice the chain has no intrinsic holding: when 9 = 1 it corresponds to 
the completely industrious random transposition walk with probability y^y to go to a 

neighboring permutation. If one inserts a holding of 1/n, i.e., P 4 - 1 + (1 — -)P, 
then the asymptotic cutoff profile stays the same and has no removable discontinuity 
at 9 = 1. 
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The Metropolis chain we defined above can be projected to conjugacy classes of S n , 
namely partitions, provided we start at the identity element. The transition matrix 
takes the following form: 



f 1- 1 A 



1A0 



(1A0- 
rfl Afi" 



V2 







if /i = A 

if /j, k = Xi + Xj 

if fa + fa = X k and fa ^ fa 

if fa + fa. = Afc and /ij = fa 

otherwise 



3. The first order phase transition at 9 = 1 for the cutoff value is not surprising, because 
the Metropolis chain has different forms for 9 < 1 and for 9 > 1. 

4. 9 denotes the inverse of the Ewens sampling parameter the chain Pg converges to. This 
choice of convention is justified by the fact that the left eigenfunctions of the chain 
Pq are the transition coefficients from the Jack polynomials with parameter 9 to the 
power sum polynomials, as derived in [12J. 

5. It will be interesting to see what happens when the 9 value in the transition probability 
is allowed to be state dependent. My conjecture is that it will always have a cutoff 
that occurs at at least |nlogn. 

The next four sections will be devoted to the proof of Theorem 12.11 



3 Preliminaries on C? mixing time 

Lemma 3.1. Given a reversible ergodic Markov chain P on a finite state space X , let fj be 
the right eigenfunctions, normalized so that 

^2fj( x ) 2 n(x) = i, 

x 

with corresponding eigenvalues 0j. Then gj(x) := fj(x)ir(x) are left eigenfunctions of P, 
with the same eigenvalues, satisfying 

*— ' tt(x) 

x v ' 

Furthermore, 

4r = E^) (!) 

71(X) 

3 

vr(x) = ;>>J(x). ( 2 ) 
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Proof. By reversibility, we have 7x{x)P{x,y) = ir(y)P(y,x). Therefore, 

ft fit?) ^p&vMv) = Y t " E r\ p ^y)^y) = J2^f\ p (y^)My)- 

* — ' z — ' TTlX) * — ' TTlX) 

y v y 

Now multiplying both sides by ir(x), we get 

Pjir(x)fj(x) = ^K{y)fj{y)P{y,x). 
y 

This proves the first part. 

The last two identities are nothing but a restatement of the fact that the matrix Q(x, y) := 
tt(x)P(x, y) is doubly stochastic, i.e., each row and column sums to 1. Here is a formal proof. 
Since {fj} forms a basis, We can decompose the function z i— > l x (z) in it: 



where the coefficients Cj are given by 
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(fj, = M^/j^M*) = fii x ) n ( x )- 



The first equality follows immediately. The second is similar. □ 

Lemma 3.2. Under the same notation as the previous lemma, one can bound the total 
variation distance to stationarity at time k starting at state x by 

4K'-^<ll^-l|||= w (3) 



7T 2 (x) 

j 

Proof. The first inequality (jHJ) follows directly from Cauchy-Schwarz inequality. To prove 
the second formula (jlj), first write 

7i ^ 7T[y) ^ 7i(y) 

Using reversibility again (in the extended form 7i(x)P k (x,y) = ir(y)P k (y,x)), we can write 

(P*(y)f _ P*(y)P y k (x) 
n(y) tt(x) 

Thus summing over y G X, we get 

P* , u2 _ P 2k (x,x) 

1 2 ~~ 7 \ l - 

7T 7l(X) 
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Next write the function y i— > P^ k (y) as the result of a row vector multiplied by a matrix: 

z 

By the previous lemma, we have 

1 x(y) = ^2c j g j (y), (5) 

j 

where Cj = £ 2 ^l x (z) gj (z) = fg. 

Finally evaluating at y — x in (E]), we obtain 

n(x) t-r 1 tx{x) 2 3 

□ 



4 Results from symmetric function theory 

First we recall from Hanlon |12] that the eigenvalues for the chain Pq projected onto conjugacy 
classes of S n , with 9 > 1, are given by 

n(A') - fl-'n(A) 

= — s — ' (6) 

For an independent proof with pointers to literature, see the proof of Theorem 14.21 

Next we derive the eigenvalues of the chain Pq, for 9 e (0, 1). Notice this is not the same 
chain as that studied in p2]. Here the identity element gets the biggest mass, whereas in 
|12j . identity has the smallest mass (Ewens sampling with parameter e (0, 1)). But the same 
result of Macdonald can be used here to derive eigenvalues. Indeed, consider the following 
matrix Tq defined by 

if a = 7i 

if a = 7f(i,j) and £(a) = £(tt) — 1 
if £(o) = £{tt) + 1 
otherwise 

where n{n) = J2 i (^) and {jii} is the partition structure of n. This quantity gives 
the number of ways to break a part in the partition structure of tc into two parts, using 
multiplication by a transposition. 

Hanlon considered the case 9 > 1 (his a is our 9), here we extend to 9 G (0, 1), which is no 
longer a Markov matrix, because the diagonal entries are no longer nonnegative. Nevertheless 
The rows still sum to 1. Then his theorem 3.5 continues to hold because the proof never 
uses 9 > 1. Likewise, theorem 3.9 holds for 9 < 1. To get Pe, we simply need to rescale Tq 



(g-l)n(Tr) 
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by 9 and add a constant multiple cl of identity matrix, c can be obtained by looking at the 
top eigenvalue. By therein 3.5 and 3.9 of [12], 9T e has eigenvalues e "( A )~"( A ) _ Thus we need 

to add 1 — 9 in order for fi( n ) to equal 1. 

Combining the two cases, we have the following formula for eigenvalues of Pq\ 

&(6>)-l-6>Al + (0vi)G) (7) 

Denote r(A) = n( - A fev^ ■ The following lemma collects a bunch of estimates about /3a: 

Lemma 4.1. Lei > 6e i/ie natural partial order on the set of partitions defined as follows: 
given two partitions represented by Ferrers diagrams X and X' , say X > X' if X can be obtained 
by successive up and right moves of blocks of X' . 

1. n(A*) is monotone andn(X) is anti-monotone in the above partial order, i.e., for X > A', 

n(A*) > n(A") 
n(A) < n(A') 

/?a is monotone with respect to the natural partial order on X. Thus fit n \ > fi x for all 
A h n, fix < /3(\i,n-Xi), and fix > fi(x u i n - x i)- 

3. Furthermore, for Xi > 

2Ai(n-Ai) 



fix<l-(9M) , , 
n[n — 1 



and in general 



/3a<1-(0A1)(1-^— 1) (9) 

?2 — 1 



In particular, if fi x (9) > 0, the above two inequalities hold with \/3\\. 
4. Finally, if fi x {9) < 0, then fi x t{9) > and \fi x \ < fix*- 
Proof, (of Lemma.) 

1. It suffices to check the first assertion for A and A' that differ by one block, i.e., A,; = 
A^ + 1, Xj = X' j -l,i< j. Then 

n(A*) - n(A") = i[A,(A, - 1) + A,(A, - 1) - A -(A - - 1) + AJ(AJ - 1)] 
= a; - Xj > 0, 

using the fact A, > Xj and A- > A'- by definition of Ferrers diagram. The antimono- 
tonicity of n(A) follows by taking transpose. 
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2. This follows directly from the previous assertion and the formula (J7|) for (3\ in terms 
of n(X) and n(A*). 

3. (JHJ) follows from A < (Xi,n — Ai) and monotonicity, after throwing away the term 

"(A) 
(»V1)(3) ' 

For we again throw away the term — in [3\ to obtain 

< 1 - (e a i)(i - ; 7 J ) = 1 - (0 A 1)(1 - ^-i). 

n(n — 1) n — 1 

4. Here we consider # > 1 and < 1 separately. When # > 1, 

Px = -^(9n(X t )-n(X)). 

If /3 A < 0, then 9n(X t ) - n(X) < 0. Since 9 > 1, n(A) > n(A*). So 

0n(A) - n(A*) > n(A) - 9n(X t ) > 0, 

which implies /3a< > |/?a| > 0. 
Next let < 1. Then we can write 



(■■) 1 (-v (■■) 

If /3a < 0, then since 9 < 1, 

n(A«) n(A) 

G3 " G) ' 



Switching A and A* we again get 



/3a* > \Px\ > 0. 



We also need some definitions and results from Diaconis and Hanlon[5]: 



□ 



Definition 1. We collect some notations to be used in the main proof below, some of which 
will be repeated. They are for the most part taken from [5]. 
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1. Given a partition Ahn, and a position s = in its Ferrers diagram (i.e., j < Aj), 
define 

h*(s) = h* x (s) := (a + l)6 + e 
h*(s) = h x ,(s) :=a9 + (£+l) 

where a = A.; — j denotes the number of positions in the same column and strictly 
below s (the arm length) , and i = A*- — i denotes the number of positions in the same 
row and strictly to the right of s (the leg length). 

2. Also define the generalization of hooklength product: 

h=h{0) :=l[K(s)h*(s). 

sex 

3. Define ca iP = C\ ;P (9) to be the change of basis coefficients from Jack polynomials to 
power sum polynomials, i.e., 

Ja(0) = 5>A, P (0)?V (10) 

phn 

4. Denote by tt = ttq the Ewens sampling measure with parameter recall 7ie(cr) = 
6~ l \°> j 'z n (6'~ 1 ), where z n (6~ l ) = Y\^ = i{0~ 1 + i — 1) is the Ewens sampling formula. Also 
let n = U n .e := ir e {l n )~ l = EKllt 1 + W ~ 

Note that when 9 = 1, j\(l) is exactly the square of the product of hook lengths of all 
positions in A, which is well known to be ( ) 2 by the hooklength formula. By Wedderburn 
structure theorem (see [9] Chapter 18 Theorem 10), we also have 

n\ = dim7r^, 

A 

therefore V\ ^tt = ^r- 

Theorem 4.2. The left eigenf unctions of the Metropolis chain Pq defined on partitions, 
normalized in C 2 (V n , 1/ire), are given by 

9xip) = (jAn/( C Xi))V2- 

with corresponding eigenvalues stated in ([71). 

Proof. We synthesize the arguments found in [T2] Definition 3.8 to Definition 3.12 and [5] 
Theorem 1. The result from [TJ5] Chapter VI, section 4 shows that the Macdonald polynomials 
are simultaneous eigenfunctions of the Macdonald operators D r q t , r = 0, . . . , n. Specializing 
to the limit q = t e , t — >■ 1 and after some afflne linear transformation, the same results 
hold for Jack polynomials and the associated Sekiguchi-Debiard operators (1431) . D$(X). The 
X 2 coefficient of this operator valued generating function turns out to be the following 
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Laplace-Beltrami type operator (our notation differs slightly frompjj] page 320): let / be a 
homogeneous polynomial of degree N in n variables, then 

D 2 e f = (-y t4 - 0V n + c n )f, 

where U n = Etxi^i? ~ *A = EIU^, V n = lE^ ^^gp 1 , and c n = 9 2 ( N 2 ) + 
6NQ) + \Q)(3n - 1) (see flS} for a proof). 

After an affine transform, we arrive at the following cleaner operator: 

L 2 e ■= Y)(\u n + l -{V n - (n - 1)N)) (12) 

which readily admits a Markov chain interpretation, when acting on power sum polynomials. 
Combining (jl5]L ( PHI) , and ( )46|) . we have 

\ As_1 

^ = ^(2^^— + )n(A) + e LyL^— ) 

Observe that for 9 > 1, the first and third terms above correspond to joining two cycles 
into one and splitting a cycle into two cycles respectively, whereas the middle term gives the 
holding probability at A. In other words, the probability of going from A to /i in one step 
under the Jack Metropolis walk is given by the coefficient of Lgp\. This translates to 

Next we show that L 2 e J\ = (3\J\, with (3 given by (jSJ) (i.e., when 6 > 1); the gen- 
eral case follows by an appropriate affine transform. In [19J page 317 it is shown that for 
the Macdonald operator- valued generating function D n (X;q,t), eigenvalues are given by 
f3\ (X; q, t) = nr=i(l + Xt n ~ l q Xi ). Now using Example 3 (c) on page 320, one can derive the 
eigenvalues for D n (X; a), by considering the limiting operator lim^i(t— l)~ n Y n D n (Y~ 1 ; q, t) 
where Y — (t — 1)X — 1. Extracting the X n ~ 2 term gives (3\, which is stated in Example 3 
(b) of page 327 as a{^)e\(a) (a is the same as 6 in our notation), since = a(^)L 2 a . 

Finally we prove the formula for the left eigenfunctions. Define the inner product (, )g by 
(PA,Pp)e = 5\^z\9^ x \ In [22] (see also Lemma 3.11 of [12J), it is shown that (J\, J\)o = j\(9) 
as defined before. Here recall the normalization of J\ is fixed by requiring that in the 
monomial symmetric function basis, its m^N coefficient be 1. Therefore expressing J\ in 
terms of p^'s, we have 

£c^ (P) =JA. 

p\-N 

On the other hand, the normalization constant for the MED(0 _1 ) distribution is z n (6 l_1 ) = 
6-\e- 1 + 1) . . . (6' 1 + n - 1) = m~ n } hence ir e (p) = Q-^fz^O)- 1 and with g x {p) given in 
pip we have 

r 2 7 r 2 e ^z 

by the previous equation. This shows g\ are indeed left eigenfunctions by Lemma 13.11 

□ 
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Corollary 4.3. The right eigenfunctions of Pq are proportional to 

Hp) = gx(p)o e(p) z p 

Proof. Since ne{p) oc 0~^^- : this follows from Lemma [3.11 and the previous theorem. □ 
Lemma 4.4. For any A h n, 

Ca,i" = !■ 

Proof. This follows from the following formula in [22], stated again in Theorem 17.11 

Ja(1 b ;0)= J] (n-(i-l)+eO--l)), 

true for all n G N, by reading coefficients of powers of n. See [5] section 4 Theorem 1. □ 

Lemma 4.5. j\ admits the following inductive bound on the parts of X: 

„• \ \ |2/i2Ai-i i^- 1 -!. -vr 2 /i26> 2 ■ 

Jx > Ai! ft A x e J(A 2 ,...,A n )- 

Note that the constant e^ 71 " 2 / 126 * 2 is not important. 
Proof. From the definition, we have 

Ai Ai-1 Ai Ai-1 

3X > [Jim II & + W(Xa,.^) = [fl(^) II (^)(! + (^) _1 )li(A 2 ,..,A„) 
i=l i=l i=l i=l 

= AJ^HAx - l)!^- 1 exp(ilogA 1 - I^)j (A2> ... iAn) 

= AJ^CAi - l)^ Al - 1 Arexp(-i^) J(A2 ,... iA „ ) (13) 

where we used the fact that 1 + x > e x ~ x2 / 2 for x > 0, applied to x = and the zeta 

sum: 

■2 



E-< 



7T 

2i 2 ~ 12' 



□ 



5 Mixing Time Upper Bound 

By Theorem 14.21 under the same notation there, four times the total variation distance of 
P^n from 7r can be bounded by 

A 
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Where we use the sloppy notation ||P* — 7r|| 2 to mean \\-^- — ly^a^. 

For A = (n), corresponding to the trivial representation on S n , and starting point 
x = (l n ), the summand exactly cancels —1: /3( n ) = 1, j( n ) = H9 n n\, and C( n vi« = 1 (by 
Lemma [4.4p . whereas vr(l n ) = II -1 , so 



7r ( 1 n)2^( 2 *)0(n)( in ) _ L 

Thus using the explicit formula for g\, we immediately have 



o2k 

lip* _ Trill = e n n\Il n V (14) 

— ' Ja 

Ahn,A^(n) 

We now break the sum according to the sign of f3\. 

o2k ± o2k o2k 

E — = E > + E — > ( 15 ) 

where J^* denotes summation skipping the top eigenvalue indexed by A = in). Next we can 
rewrite the first summand on the right according to the size of Ai, and obtain the following 
bound: 

o2k 1 o2k 

e n n\u n y tL= y y «nA 

s=n-l A:Ai=s l( 3 A >0 JA 

< £ max{/3f : & > 0, A x = ,} V (16) 
— — ?x 

s=n-l A:Ai=s JA 



Splitting IT = Ii n e into two subproducts, and using Lemma 14.51 we have 

Tln\9 n < - n": Al (n-l/ " 1 n w _ A > - AQIfl^ 
Ja " ^^AiPA?" 1 "^-^/^ 2 J(a 2 ,...,a„) 

where the second quotient factor happens to be ^J^^ ff^ a„)(1™~ Ai ) ( see Theorem 14.21 
and Lemma EOJ. Also denote the first factor by g n ,Ai- 

By d2]), the definition of IT n _ s := ^(l™- 5 )- 1 (see' Definition [TJ, and the fact U n = 6 n (n - 
l)\e e 2 -' i = 2 ~, we can bound the summand of the RHS of f|T6|) for a fixed s as 

' 7\ ' 7T l n s r p 7r 1" s s si [n — sy. 

We will now reduce the C. 2 bound (1141) to bounding the following quantity: 

k + ■= e ^ + i-r ^ n - s ^fj — v' ^ 

z — ' ' s s! (n — s)! 

s=n— 1 v ' 



12 



where f3 s>+ := max{(3\ : (3\ > 0, Ai = s}. 

For the second summand of (1151) . we obtain 

J2^<t ■ & < °' M = 4E + (is) 

/3a<0 JA «=n-l A:A*= S 3X 3lU 

Using the explicit formula (jTJ) for /3a, we get 

f3 ln = 1 - (0 A 1) - (r 1 A 1) G (-1, 1), 
for ^ 1. On the other hand, by Lemma [4.51 

n— 1 n— 1 Ti—l n— 1 

UnnlP/jm = + i0)n!0"/(n! + *)) = * + II^ + *)> 

8=1 j=l i=l i=l 

which is of polynomial order in n for each fixed 9. Thus 

0^ n) n n nW n /j ln = o(l), 

which has negligible contribution in f fl4|) . 

For the remaining terms in ( fl8l) . first observe that, 

Ai Ai-l 

i=l i=l 

hence when 9 < 1, 

Ai Ai-l 

3X< > U( i9 ) II ( W + ^'(W-)* 
i=l i=l 

and using the bound |/3a| < /3a*, for (3\ < 0, we get 

3x Ja 
If > 1, the j\t is comparable to j\ within an exponential factor: 

jx* > Jx9 2n . 

Furthermore by the explicit formula of /3a, we have 

o (fl s M^) ,- X n(A) x n(A) n(A') x n(A)j-n(A^) fl . lfl m 

So since /c = 0(nlogn), Q- 2k Q 2n = o(l). Thus we can still compare the negative /3a sum to 
the positive one: 

I3 X <0 JX 
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It remains to bound 2b n+ . First note that the factor 2 in front is immaterial, since for 
A ^ (n), 

0X < /3(n-l,l) < e_Qe(C) > 

thanks to the monotonicity of flx's. So by increasing c in k = 2 (6>/\i) n ( c l°& n )' we can 
decrease 6 nj+ by a factor of 2. The factor eg = e 7 ^/ 26 * 2 can be ignored similarly. We can also 
get rid of the factor _1 in ffTT|) as follows: 

For a > n/2, (H)^ 1 -! = 0e {\), so again increasing c annihilates it. For s < n/2, recall 
the second bound on f3\ Qj, which implies that for Ai < n/2, f3\ is bounded away from 1 
uniformly in n. Now in the definition of b n>+ , (3\ is assumed to be nonnegative (alternatively, 
0\n is bounded uniformly away from —1), hence raising (3\ to the power Q(n) easily cancels 
any power of n, i.e., 

r r l - i p? = o{i). 

So together, by increasing c, we can reduce the problem to bounding the following quantity 

i ns n\ 

s=n—l 

The only estimates we rely on now are (jSJ) and (Q from Lemma I4.lt the idea will be 
similar to [8] (see also [7]). First note that it suffices to show 

Qn—s | 

^) 2 /3f + = 0(1), (19) 



[n — sjl s 

uniformly for all s6 [l,n — 1] and c sufficiently large. Indeed using ([9]), we have 

S ,+ < e~ x ~ x2/2 < e~ x 

for x = (9 A l) 22 ^. Therefore 

1 n 

0(1) ffi <o(i)j2^ tc = o c (i), 

s=n— 1 t=l 

by geometric summation; in fact, using (jSJ) we can get a better bound, but that's not 
necessary. 

Next recall (J5]) as well as the estimates (no Stirling formula needed) 

^" ^ J™ log xdx+log n— log s n log n— s log s— (n— s)+log n— log s 

s! ~~ 

( n _ s )| > e Sr~ aio S xdx = e (n-s)log(n-s)-(n-s-l) ^ 

Taking logarithm, we can bound the left-hand side of ffl9l) by 

log[- rr(-r) 2 /3s,+] < (1 _ 2a)(l - a)nlogn - (1 - a)nlog(l - a) - 2cm log a - 2 log a 

+ (Ci(0)-2caO(l-a)n + C r 2j (20) 
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where C\{6) is a constant that depends only on 9 and Ci is another constant. 
For a > a G (1/2, 1), the right hand side of (12"U|) can be further simplified to 

(1 - a)n[{\ - 2a) hgn - log(l - a) + C[{6) - c] + C 2 . 

Since (1 — a)n > 1, and we can choose c as large as we want, it suffices to show the expression 
inside the square brackets above is 0(1). But when a = 1 — l/norl/2, this is clearly true. 
Furthermore, the derivative 

— [(1 — 2a) hgn — log(l — a)] = —2 \ogn H — 

da 1 — a 

is monotone increasing, showing that the function a i— y (1 — 2a) hgn — log(l — a) is convex 
and its value for any a G [1 — 1/n, 1/2] is bounded by the values at the boundary points. 
Next let a < a®. Using the second bound for {3\, ([9]), we have 

Then the logarithm of the left-hand side of ffl9|) is bounded by 

(1 — a)n log n + (1 — a) (log 9 — log(l — a) + l)n — 2an log a 
+21og« - (1 - a)n\ogn - (1 - a)cn < (1 - a)(C(6) - c)n + 2(l - an)loga. (21) 

Clearly (1 — an) logn = 0{n) for a G [-, oto]. So for sufficiently large c, the right-hand side 
above goes to — oo. Together this shows f|T9|) is true for all s G [-, 1 — -], and concludes the 
upper bound for the mixing time. 

6 Mixing Time Lower Bound 

We rely heavily on results from [22] . Again we collect some notations needed in the analysis 
below: 

Definition 2. For X,p\- n, let 

• H\ be the product of all hook- lengths of the Ferrers diagram for A, 

• z p :— niLi i mim i- and nii = nii{p) is the number of parts in p of length i. 

• Xx(p) De the character of S n indexed by A evaluated at an element of conjugacy class 
p; alternatively, they can be defined by the system 

Pp = ^2xx{p)s x , 
x 

where s p are the Schur polynomials. 

Warning: Note the Schur polynomials are not direct specializations of the Jack polyno- 
mials; they differ by a factor: 

SA = # A -V A (1). (22) 
Thus the matrix Xx(p) is inverse of H^cx^l). 
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Lemma 6.1. ([19] Ch. I, (7.5); see also [22] (50)) The relation between c\ tP (l) and Xx(p) 
is given by 




p 




p 



Comparing coefficients with ([TO]) in Definition [JJ yields the result. 



□ 



As in the 8 = 1 case studied by Diaconis and Shahshahani, the strategy will be to use 
a certain eigenfunction / of the chain as test function and compare the probabilities of the 
event {/ < rj} for some i)6K under the stationary distribution and the distribution at time 
slightly before the mixing time, which in our case is k(c) := ^W n(logn — c). 

In the case where 9 = 1, p i-> X(n-i,i)(p) — m i(p) — 1 is the desired eigenfunction. So 
it is natural to guess that a suitable affine transformation of the fixed-point (aka 1-cycle) 
counting function p i— > mi(p) is the desired eigenfunction. 

Lemma 16.11 shows that the following normalized version of c\ tP is the right analogue of 
characters of the symmetric group 



Thus our candidate test function will be d\(p) = d\ tP . 

It is straightforward to compute Poo(^a < v) where denotes the stationary measure; 
the number of cycles are asymptotically independent and Poisson distributed. To estimate 
Wk{d\ < i]), one uses the second moment method. The first moment of d( n -i,i) is easily com- 
puted since it is proportional to the right eigenfunctions of the chain Pg (see Corollary 14. 3p . 
For second moments, we need to decompose c?L_ 1 ^ as linear combinations of other d\s. 
this is accomplished by first expressing d( n _i,i) and other c^'s in terms of powers of m^'s, the 
number of i-cycles (not to be confused with monomial symmetric functions), then deducing 
the relationship by solving the appropriate system of linear equations. The analysis below 
will be an elaboration of this strategy. 

First we need 

Proposition 6.3. ([22] Proposition 7.5) 

Let d\ )P = d\ tP (9) be defined as above. Then 



dx tP {9) := c x , p {9)z p 9-^-^H^. 



(23) 



fc 



6 k+ \k + n)d (n ^ p {6) = ^(-l)^(j + („ + k 



3=0 
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Note that the partitions in the proposition above is for n + k, rather than n. 
From this, we easily obtain 

i-i,D,^) = ~\ + 1 + ( g~ 1)g ™i(p), ( 24 ) 

and 

, n \ 1 ,1 + in -1)9 2 + [n -2)9. . . 

2 + (n-2)fl , .„ 2+(n-2)# . . 



Note in particular, 



X(n-i,i)G°) =mi(p) - 1 (26) 
3 1 

X(n- 2 ,i 2 )(p) = 1 - 2™i(p) + 2 m i(p) 2 - mz(p), (27) 



as expected. 

Using the Schur-Weyl relation 



we also obtain 



Xn-IA — Xn + Xn-1,1 + Xn-2,2 + Xn-2,1 2 , 



1 2 3 

X(n-2,2) = - 7j m l + m 2- 



To get J( n _2,2)) w e need the conjecture right after Proposition 7.2 as well as Corollary 
3.5 from [22] • The conjecture has been proved in [IB]. Notice the parameter a is the same 
as our parameter 9. 

Proposition 6.4. ([22] Proposition 7.2) The Jack polynomials corresponding to the partition 
(2 l , V) has the following expansion in terms of the monomial symmetric basis m\: 

i 

J(2\v) = ^2{i)r{9 + i + j)r(2(? - r) + i)!m(2r jia (i-r)+j), 

r=0 

where (i) r :— i(i — 1) . . . (i — r + 1). 

Proposition 6.5. (ffflf Theorem 1.1) In terms of Schur polynomial basis, we have 

% 

J(2\v) = 22(i)r(9 + i+j) r (i-r- 9)i- r (i + j - r)!s (2 r ) i2(i-r)+j). 

r=0 

The next result relates Jack polynomials corresponding to conjugate partitions, when 
expressed in terms of the power sum basis. 
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Proposition 6.6. ([22] Corollary 3.5) Let J x = ^ c\^{0)p^, then 

We also recall from Corollary 16.21 that s\ = Y^ p X\(p) z p 1 Ppj where Xa(p) is the character 
of A evaluated at p. 

Combining the previous results, we easily get 

J n -2M = 5>0r^[(2 - - i)(„ - 2) X J" +2(n-2+ 1 -) 
p 

(1 - i)(n - 3)!x2i,i»-'(p) + 2(n - 2 + i)(n - 3 + i)(n - ^^(p)]^, 

where X\{p) is the irreducible character A evaluated at p. Therefore we can read off the 
coefficients: 

<VW*) = K —^ p ((2 - f )(1 - tX™ - 2)! X m(p) 

+ 2(ri - 2 + i)(l - - 3)! X 2M»- 2 (p) + 2(n - 2 + - 3 + ^)(n - 4)!x 2a>1 »-«(p)), 

(28) 

using the relation Xx t (p) = XA(p)sgn(p), as well as the formula for Xn, Xn-i,i, Xn-2,2 and 
Xn-2,1 2 derived above, we also get 

n — 2 + - 13 1 

VW') = ">l (w _ 1)(w _%) [(" " 3)(1 " ? ) - ^ " 3 + ? )] 

, , (n-2 + i)(n-3 + i) , m (n-2 + l)(n-3 + i) 
+ mi 2(n-l)(n-2) + ™ 2 (n-l)(n-2) ' (29) 

Using i = and j — n, one gets 
hence 

W#) = Q n ~^n\z-\ (30) 

and 

d ( „),p(0) = 1. (31) 

To mimic the case of 9 — 1, we need to express d( n -i,i),p(#) 2 in terms of the other o? A 's. 
First using ( 12~4|) . we get 

l On n On n 
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Using mi,m^,m2 and 1 basis, we can write 

d(n-i,i),p( ) 2 = u + ud( n _i j i) >p (0) + «;d( n _2,ia) iP (0) + xd {n _ 2 ^) )P {9) 2 , (32) 



for some indeterminates u,v,w } x. 



Comparing coefficients of the mf s, we get the following four equations: 

^ + ^ = ^ ( 33 ) 
.1 77-1. ,l + (n-l)0 2 + (77 -2)0. 
"fe + ~ } " " ( ^ + 29^- ) 

77-2 + i r/ w 1 3, L. 2 1 77—1. 

+g (n-l)(n-2) K W - 3 ^ 1 "^-2^- 3 + ^ = S^n + ~ > ^ 
2 + (n-2)fl (77-2 + i)(77-3 + |) _ 1 77 - 1 
W 20n +X 2(n-l)(n-2) ~ [ 9n + n j 1 ' 

_ w 2_±(n-2)9 (n-2 + l)(n-3 + i) _ 
W 0% +X (n-l)(n-2) " ( j 

Solving, we get 

_ (77 2 - 477 + 3)0 3 + (n-l)0 2 + (ri-l) 2 + 77 -3 _ 1 1 1 

W ~ 2 (0 + 1)(0(t7-3) + 1)t7 ~ 2 (0 + l) + 0+T + n 

(t7 3 - 6n 2 + lira - 6)0 3 + 2(2r7 2 - 7n + 4)0 2 + (3r7 + 2)0 - 4 _ 1 
W ~ 0n((n 2 - 5n + 6)0 2 + (3r7 - 8)9 + 2) ~ 1+ V 

2(l + 0(r7-l)) 2 20 1 

U ' (1 +0)(2 + 0(t7-2))t7 1 + V 

2(l + g(77-l)) 2 (n-l)(r7-2) 2 I 

(l + 0)n 2 (l + 0(277 - 5) + 2 (n-2)(77-3)) 1 + V' 

Notice x + w = 2 + 0(£). 

Also by (ED, (P25D, and (ED, we get 



dn,l n — 1 

(ra - l)(ra - 2) 



^(n-l,l),l n — n—1 



(n-2,l 2 ),l™ 2 

77(77 — 3) 

«(ra-2,2),l™ = ^ • 

which are independent of 9, because of the normalization chosen for Jack polynomials. 
Finally we recall c\{9) are eigenfunctions of Pq, hence so are d\(9), with eigenvalue (3\. 



19 



We list the relevant eigenvalues here: 



0(n) 


= 1 




£(n-l,l) 


= 1 - 


(0A1) + 


?(n-2,l 2 ) 


= 1 - 


(9A1) + 


P(n-2,2) 


= 1 - 


(9A1) + 



nln-2\ 2 4 

l-(0Al)- + (9(-). 

n n z 



2 



(fvi)© 



Notice 



hm ( 1 -fr-l.l)) = lim (l-^(n-l,l)) =2- 
ra-S-oc 1 - ^( n _ 2j i2) n->oo 1 - /3( n _ 2j 2) 

Now we are fully equipped to prove the lower bound. First observe >Cx>(^(n-i,i)(0)) ~< 
Poi(0 _1 ) + 1, which comes from Feller coupling. Here d( n -i,i) stands for the random variable 
d(n-i,i),p where p has Ewens sampling distribution with parameter 6*" 1 , as indicated by the 
subscript oo. therefore 

lim P oo (d (n _i 1 i)(0)<77) = l- 

r/— >oo 

Furthermore, 

vax fc (tZ (n _i ) i)) 



P fc (d(„-i,i)(0) <»7) < 



(?7 - Efe^!,!))) 1 



Let k = \(9 1 V l)n(logn — c) for any c > 0. 

Since g?a are eigenfunctions, we can compute the mean and variance at time k: 



Efcd(n-i,i) = (n - 1)(1 - (0 A l)-) fe + 0(1) = e c + 0(1) 



var fc d (n _ 1;1) = E fc 4,_ lfl) - (E fc rf ( „_ l5l) ) 2 = « + (n - i^-^Df 

+ ^ (w- l)(w-2) w + w(w-3) ^ e _ (gA1) ^ (i+o(i)) _ ^ _ ^2 e -(eAi)f 

Therefore if we let rj = |e c , then 

e c 

lim . limjnf P fe ( c )[d(„-i,i) < rj] < lim lim inf 0(1) - ^ - = °" 



c— >oo n c— >oo n 



Thus 



lim lim ||£i»P fc(c) - 7t||tv > lim liminf |F fe(c )[d ( „_i,i) < ^e c ] - Poo[d(n-i,i) < ^e c ]| = 1. 
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Remark 2. 

Wilson's method gives a suboptimal lower bound, once we know the "geometric" information 
that d(n-i,i), P = _ \ + 1+ ^on 1 ^ m i(p) : Let X x be the random variable distributed as 5 X P. We 
have 

R : = sup E(d (n _ 1>1 ) jXl - rf (n _ lil)j:c ) 2 < 1, 
ie5n 



log-J— = (0A1)- + O(^ 



^(n-1,1) 



n n 



2-" 



and d^-i^j^n — n — 1. Hence by Wilson [18J, 

< - hW - 2 1 og(i/A„-»,) [ '° sl m 1 + log — 1 

> 4(1 A g) logn + loge^ + OCl). 

This misses a factor of 2 from the lower bound obtained by second moment method. The 
discrepancy is possibly due to the non-local nature of the random transposition walk. 



7 Metropolized random transposition walk starting from 
the n-cycle class 

In order to use the machinery of Jack polynomials to diagonalize our Metropolis chain, the 
starting state must be a union of conjugacy classes of S n , i.e., the chain we are considering 
effectively lives on V n . 

Stanley proved the following conjecture of Macdonald and Stanley: 

Theorem 7.1. (JMj Theorem 5.4) 

Consider Jack polynomials with parameter 9 (suppressed), J\(x) where x G M™, A h n, 
restricted to the first n arguments, under the normalization P(3) of Theorem 1.1 in [22], 
which is equivalent to the transition coefficient c\^n = 1 for all X, then 

Jx(i n )= n (n-(i-i) + e(j-i)), 

where the product is taken over all blocks in the Ferrers diagram of X, specified by the usual 
coordinate 

As pointed out in [22] and [5], this formula allows one to recover the coefficients C\ jfJ , for 
some special /x, for which the number of blocks £(n) uniquely determines /i. For instance, 
/i = (1™), (2, l n ~ 2 ), and (n). For generic n, the list above is exhaustive. 
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Note that p^(l n ) = n e ^\ which is always a monic monomial in n. Thus we have 

c Aj i« = 1, (37) 

C A,(2,1™~ 2 ) 



CA,(n) 



E - 

(M)eA 

n 

(i,i)#(i,j)eA 



r 



(i-1) = 0n(A*)-n(A), 
(i-1)]. 



(38) 
(39) 



Recall that the left eigenfunctions indexed by A evaluated at fi is given by 

C A,M 



9M 



(i A n/(^n!))V2' 



where II = 9"( 



j. /i will be the starting point of the chain. We first exploit c\ / 2 i 



Recall the formula Ti(fi) = 9 



-1 n! 



Thus tt((2, l n - 2 )) = fl^Sfef. Using the upper 



bound lemma with left eigenfunctions, and the fact that 7r(2, 1" 2 ) 

,dP fc 



^©/(On, we get 



4 P K 



7T 



TV 



< 



1 1I 2 



tt(2,1™- 2 ) 2 



<7a(2, 1 



n-2^2 p2k 



[0 



l-n dn ( n - !)i 



E 



[0n(A* 



jAn/(n!^ 



<A)] 2 
"Pa 



1. 



(40) 



When # > 1, the stationary measure puts more weight on (n) cycles and less on the identity 
element. So the Metropolis walk moves from (l n ) to (2, l n_2 ) with probability 1 in the first 
step, and we expect the mixing time to be essentially the same as if the walk starts at l n . 
This is indeed the case as we demonstrate below. The A term of the last expression is 



fln(A') -ra(A) 



n\9 n U 

Jx 



fix 



2 k 



When A 



n 



1,1), we get (1 - Oin- 1 )) 



nW n U o2k 



Pi. 



which matches the leading term 



'J(„-L1)>-1,1)' 

of the corresponding bound for the walk starting at 1™. Next observe that for 9 > 1, 
|0n(A*) — n{X)\ < 1 for all A. Thus each summand in ( 1401) is dominated by the corresponding 
term in the upper bound formula for the chain starting at (l n ) (fl4|) . Since each summand 
in (111)1 is positive, we get the same mixing time upper bound for P(2,i«-2)- 

When 9 < 1, we see that as long as 9n(X t ) — n(A) > 0, the terms in (I40p are still 
dominated by the corresponding terms in (fl4l) . For A with 9n(X t ) — n(X) < 0, \i3\\ is 
bounded by (1 - 9 A 1) V |1 - 9 A 1 - 9' 1 A 1| = 1 - £1(1) provided 9^1. On the other hand, 



the factor 



6>n(A i )-n(A) 



is at most 9 2 , which can be easily compensated by increasing c in 



the mixing time t(c) = gj^n^c + logn). 

More interesting is P( n ), the Metropolis chain starting from the uniform measure on all 
n-cycles. 
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Theorem 7.2. The mixing time for P( n ), T m ix(|) < @e{ n ); f or 9 — 1> an <^ ^ e version 
of the chain Q = ^^-P + we nave £/ie following more refined upper bound 

. . n , . 1 — 4e 2 . 
r mix (e)< I log(^^). 

For the lower bound, we only consider the case 9 = 1: 

Theorem 7.3. For = 1 and any sequence t n = o(n), 

\\P^-7r e \\TV=l-o n (l). 

Remark 3. The lower bound does not match the upper bound exactly; getting that pre- 
sumably requires higher moments of the test function involved. However it does show that 
Q(n) steps are necessary for mixing. 

Before proving the above result, we need a definition and a lemma. 

Definition 3. The 9-Durfee height of a Ferrer's diagram A is given by min{w > 1 : 
(u, \u9~ 1 \ + 1) ^ A}. Intuitively, u x (|_0 _1, wj + 1) is the maximum size of a recangle 
with aspect ratio 9 that could fit into A. This is a generalization of Durfee square, which 
corresponds to 9 = 1. 

We also call the u x ( |_0 _1 wj + 1) rectangle in A its 0-Durfee rectangle. 

The following lemma is crucial in the proof of the theorem; I thank Bob Hough for 
providing one of the key ingredients. 



Lemma 7.4. Let s = \\. Then 



c 



2 

A,(n) 



e 



[Ci(e)(n-s) 1 /2+C72(e)] log(n-s+l) 



h ~ C,{9)s^ ' (41) 

for d(0) = 9 1 ' 2 + 9- 1 ' 2 , C 2 (9) = 1 + 0- 1 and C 3 (6) = 9e~^^ 12e \ 

Proof. The first observation is that we can divide A into two parts, 
Consider the following three sets: 

A+:= {j eN:(i,j)e\,F(i,j)>0}; 
Aj:={ieN:(i,j)eX,F(i,j)<0}; 

A° :={(i,j)e\,G(iJ) = 0} 

where F(i,j) = [9(a(i, j) + l)+£(i, j)][9a(i, j)+£(i, + is the (i,j)th factor in the definition 
of j\. Also define G(i,j) = [9(j — 1) — (i — l)] 2 , the factors in c\ ,y 

It is clear that for every G A, there are three possibilities: i G A+, j G A~, or 

(i,j) G A . The last one includes the block (1,1), which is excluded from the defining 
product for C\^ n y Let's also define A+ := minA+ and \J := minAj. Then we can match 
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off F(i,j) and G(i,j) according to the partition above. First we consider the first row of A, 
which consists of Af, for both G and F, together with F(l, 1) in the denominator, 

f q-> n wfT) = n (e(j -v-v- i ^ 2 ' n - n + vw* - * + ^ 

J>2 ^' ' j=2 j=l 

s-1 

= 02(»-i)( a - l)! 2 /(^- 1 s!( s - 1)! JJ(1 + fl-^- 1 ) 

i=i 

1 



ese 0-l[log( s -l)-2j] 

This provides the first factor in the bound (j4"Tl) . 

We can now partition the remaining blocks in A into the above three sets. Consider first 
Af for i G [2, h], where h is the #-Durfee height of A: 

II J^T) = ft 00" - !) - (* - ilf/HW-i) + !) + + '(U) + 1) 

A.-A+ + 1 A,-A++l 

< n w)v n 

3=1 3=1 

< ^ |A?l (|A+l07^ 2|A?l |A+|!(iA+| - 1)!] = |A+|. 



Similarly, we have for j > [1, [hd x \ +1], 

n^<iA7|. 

11 F(i,j) ~ 1 > 1 

Next observe that |A+| < n — s for i ^ 1, and |Aj | < n — s + 1 for all j, since besides the first 
row, which has length s, the other rows and columns all have lengths bounded by n — s + 1. 
Summing up, we have 

c 2 , , h LW _1 J+i / 

f^< (to «-')-niAri n w^^^-, — ■ 

Jx i= 2 i=i ys 

where /i is the #-Durfee height of A. To bound h, simply observe that (h — 1)( [_(h — 1)# _1 J + 
1) < n — s, since the biggest possible Durfee rectangle is achieved when the rows excluding the 
first one form exactly a rectangle with aspect ratio 9. From this we have (h — l) 2 ^ 1 < n — s, 
hence 

h + he- 1 < V^^(o 1/2 + e- 1 ' 2 ) + i + e-\ 

This accounts for the second factor in (j4"Ti) . □ 

We also need to bound c\ i n Jj\ in terms of the transposed tableau A T , to account for the 
negative eigenvalues. This is provided by the following corollary: 
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Corollary 7.5. For A with \{ = s, we have 

where C\{9) and C2(9) are defined as in the previous lemma, and C' 3 (9) = Oe^^ 9 ^ 12 . 

Proof. The numerator in the bound follows in the same way as the previous lemma. For the 
denominator, which comes from matching F(i, 1) for i G [l,s] and G(i, 1) for i G [2,s], the 

product Yli=2^(^j) = r3i=2 ( — (* ~~ -0) 2 = ( s ~~ l)' 2 doesn't depend on 6* and the product 
rii=i F{h 1) can be bounded below by 

s 

Y[(s - i + l)(s -t + 6)> s\(s - i)!0 e 0i°g(*-i)-T 2 <Vi2_ 

□ 



i=l 



Proof, (of Theorem 17.2 j) 

As in the case of starting the chain at 1™, we have the upper bound formula: 

^WPU ~ Alv < \\Pl) ~ Al 



^ y- c l,(n) n\9 n 2t 



A 



First we look at the case 9=1, where we insert a 1/n laziness, i.e., use the following 
kernel ^^P + -I- In this case, c\ t ( n ) = for all A except the ones that look like a hook. 
That is because as long as the block (2, 2) is present in A, the product defining cx,( n ) will be 



c 2 



zero. For A = (Ai, l n Al ), we have A ^ n) = n 2 . This exactly cancels ^77-^2 • Furthermore, 
the eigenvalue is given by 

_ n- 1 n(A t ) - n(A) 1 

PA 7n\ r -• 

n g) n 

/V\ /n-s + l\ 

Notice that /3 (Sj i«- s) = n ~ U2; ^ 2 ; +^, from which it is easy to see that for s G [n/2, n— 1], 

> |p\n-s+i,i 8_1 )l- Thus modifying the bound in Lemma |4~T1 to account for laziness, 
we have for Ai > n/2, f3\ < 1 — 2Al (^~ Al ) _ Therefore for t = cn, 

1 



^ - 1 (2) -J_2 ) , 1 l2t t f 1 - 
I") 

s=n-l \2J 



< 2 v [ y2 > " 2 ' + -r + (-) 

n (") n n 
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n/2 

E n _ l _ 2s(n-s) 1 
n 

s=n— 1 

< 2e- 4c /(l - e" 4c ) + o(l) + e 
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Thus the mixing time is of order 0{n). In fact by introducing more laziness, say of size 
log n/n, we can eliminate the factor of 2 and the (l n ) term e~ c in the last expression, and 
achieve a total variation upper bound of e _4c / (1 — e _4c ). So we have to solve 

e ~/(l - e~~) < Ae 2 . 



This gives r mix (e) = f log(^-,. 

Similar argument can be used when 9 is a rational number, say p/q. In that case, any 
A h n with (p, g) e A has cx,( n ) = 0, hence by considering the 6*-Durfee rectangle, one has at 
most s p+q A's with n — X\ = s whose cw n ) is nonzero. For each such A, a weaker version of 
the argument in Lemma EH gives that c\ t n J j\ < s p+9 . So one has to bound 



4st 



£ 

s=n— 1 

which is < 1/4 when t > + q). Thus the bound grows linearly in the degree of irra- 
tionality of 6, using this argument. 

For general 8, one needs the full strength of Lemma 17.41 and Corollary 17.51 Also recall 



n!0" _ (9-%n* < 



vr((n)) 2 n nlfl" 2 

Then using the Hardy-Ramanujan asymptotics for the number of partitions of n , \V n \ < 
C(e n ), as well as Lemma |4TT| the upper bound formula becomes 



I*-) -II' £ A £— + JE JC — + — 



i 

1/2 OCfn-.^V2l loB-Cn-^n - 2a ^"— j* x , , /.-lN2t 



< 0n 1+ ^(2 ^ 0(l) e ("- s ) 1/2 e°«"- s ) 1/2 ) lo ^ s+1 )e'^t) + (61 A 

s=n— 1 



s=n— 1 

Thus setting t = cn, the last expression clearly goes to as c —¥ oo. 

□ 

Next we deal with the lower bound. The main difficulty is that if one takes the number 
of fixed points, m 1 = 1 — x™ -1 ' 1 ; as the test function, the best one can get using second 
moment method is 

P^X"" 1 ' 1 = 1) - PooO^ 1 ' 1 = 1) > 1 - e" 1 + o n (l). 

This is essentially because the starting value of m 1 at the n-cycle class is too close to the 
equilibrium value. So in order to distinguish <5( n )P <n and ng better, we will use the following 
sequence of test functions m\ +rri2 + ■ ■ ■ + rrik, with k going to oo. First we need to relate 
the mj's with irreducible characters x X °f S n - 



26 



Proposition 7.6. For all k < n, and v h k with part lengths all < n — k, and let (n — k, v) 

he the partition consisting of the union of (n — k) with v, we have 

• \ n v is a linear combination of cycle counting functions (not to be confused with 
monomial symmetric polynomials) m M := YYi=2 m m> w here p are partitions of n, with 
part lengths arranged in nonincreasing order and pi>n — k (note we start at i = 2). 

• Such representation of x n ~ k ' v is unique. 

• {m^ : fii > n — k} is an independent set spanning {x^ '■ Pi > n — k}; in particular, 
{mi, . . . ,rrik} is independent for k < |. 

Proof. We record the following explicit formula for characters of ^([TH] Chapter I §7.8): 

X X p = [x X+5 }(a sPp ). (42) 

The nice thing about this formula is that when A is close to (n), such as (n — k,l h ) for 
small k, the number of terms 0,5 = ^2 weS e(w)x wS that contribute to the coefficient above is 
small, hence the right hand side can be expressed as some low degree polynomial function 
of mi(p), m 2 (p), . . . , m k (p) where m^p) denotes the number of parts of length i in p. Using 
(|42p . it's easy to see that x n ~ k ' u is a linear combination of terms of the form [xj 1 . . .x] a ]p p , 
where Y2i=i 7« = n i Ti — 71 ~ ^ anc ^ ^ < 1 + E(v). Thus it suffices to show such a term has 
the required representation in terms of the m M 's. 

Let A 2 h 72,..., X s h 7 S be a sequence of partitions. Each such sequence leads to 
one way of extracting coefficient from p p , by matching each element in the set 

{x Xj : 2 < i < s,j G [^(A*)]} with a distinct factor in the product p p = rii=i£W So we can 
further reduce to showing that the resulting contribution, denoted [71, A 2 , . . . , X s ]p p , has the 
required representation. 

First we combine the A 2 ,...,A S into a big partition A h ^2 s i=2 li < k. Write it in 
exponential notation as A = (l ai , . . . ,t at ), i.e., t is the longest part length in A; in particular 
t < k. Furthermore we write aj = ^2j=i a ij so that ^ = • • • > t ati )-, an d let 

m(p) \ ( m(p)\ ( mi(p) - acn\ ( mi(p) - Yl%\ a ij 



oca,...,oc it J \ ocu / \ a i2 J \ a it 

be the multinomial coefficient of rrii(p) choosing the a^s. Then the contribution of the 
xj 1 . . . x] s coefficient in p p coming from the sequence A 1 , . . . , A s is given by 



[ 7l ,A 2 ,...,A s ]p p 



m i(p) \ ( m i(p) \ ( m t(p) 
a n , . . . , a is J W, ...,a 2s J \&ti, • • • , oi tl 



Notice that this is a polynomial in mi(p), . . . ,m t (p) such that for each term of the form 
mf 1 . . . mf', Yfi=iifti — So it can be written as a linear combination of cycle counting 
functions m p with pi > n — k. 

Since x M are orthogonal functions on S n with respect to the uniform distribution, they 
form an independent set. This immediately implies the independence of {m p : p\ > n — k} 
by cardinality consideration. The uniqueness of the representation of \^ a l so follows from 
independence. Finally by definition of m^, m k = vn n _ k ^. □ 
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The following result might be of independent interest: 

Corollary 7.7. Let C(S n ) be the space of class functions on S n , and let ch : C(S n ) — > A n 
be the characteristic map given by x X ^ s a? the Schur polynomial indexed by A. Then the 
image of {m^ : /x h n} under ch gives a basis of A n , the space of symmetric functions of 
degree n. 

Proof, (of Theorem 17.3 j) 

Let f k = m 1 + ... + m k and write f k = Eai-™ c aX A - By ortho normality of x A , we have 

E c a = Wf k ~ c (n)ll£2(7r) = var^(/fc) = log k + 0(1), 

(n)^Ahn 

where the last equality follows from the strong asymptotic independence of mi, ... , m k under 
the uniform measure 7r and that m.j converges weakly to a Poisson variable with parameter 
i^ 1 (see pQ). Furthermore we have x\ n ) 7^ if and only if A is a hook, that is, A = (n — k, l k ) 
for some k, in which case x\ n ) = this follows for instance from Lemma [6. II and (|39|) . Now 
we can estimate the mean and variance of f k at t = t n : 

Etn(/fc) = E C ^AX(n) = E C **(») + E - l)cAX(n)- 

A^(n) A^n A^(n) 

The first term above equals mi + . . . + m k — C(„) evaluated on an n-cycle, hence equals 
— c (n) = — Eoo( m i + • • = — Ei=i The remaining term can be bounded in absolute 

value by 

E icAX? n) i(i - = E [<=aI(i - ^) < E M 2(n ~ Al) * (i + o(i)), 

A A A 

where we used the estimate f3\ = 1 — 2 ( n ~ Al ) + (9(n~ 2 ) (see (JTj) with 6* = 1) and the fact that 
t = o[n) so that for fixed k, 2(n ~ Al ) t = o(l). By Cauchy-Schwarz inequality, we have 

c A — < E N ^ v^OogA;) 1 / 2 = o 1 , 

n L — ' n n 

A A 

which shows 1E tn (f k ) < — log k + o(l). This constrasts with E OG (/fe) = 0. It remains to bound 
the variance var in (/ fc ). Let /| = J2x\-n c 'xX X f° r some c' A . Then c' x = unless \\ > n — 2k; 
this can be seen by expressing x X m terms of m M . Using x A n ) £ {0, 1} again, we have the 
following estimate 

Var t / fc = E^Wn) - ( E C ^iX(n)) 2 = E C ^W ~ ( E + R * = ^> 

A A^(n) A A^(n) 

where ^ = £a c'M - l) X A n) - El c xX \ n) El c A (# - l) X \ n) ~ (El c A (^ - l)xf n) ) 2 . Since 
I El c aX A „) I < l°g fej we get -Rfc = o(l) by the same reasoning as before. This shows \ai tn f k = 
o(l) and f k is essentially constant at t = t n . By Chebyshev inequality, we have 

V tn (fk < -c (n) /2) > 1 - - = 1 - o(l) 

(-c (n) /2- (-c (n) ) + o(l)) 2 
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On other other hand P oc (/ fc < -c (n) /2) < < + o(l) = 4/logfc + o(l). 

Therefore 

liminf ||P ( *;\ - vt||tv > liminf liminf P tn (/ fe < -c (n) /2) - Poo(/fc < -c (n) /2) = 1 - o(l). 

□ 

Remark 4. Theorem 17.31 can be easily extended to 6 ^ 1 by replacing x~ with d\~. The 
analogous result to Proposition I7.6l can be derived by precomposing the transition coefficients 
from x X to m M with the transition coefficients from Jack polynomials to Schur polynomials. 
The latter form an unipotent (upper triangular with ones on the diagonal) matrix under any 
total order compatible with the natural partial order on V n , because it can be expressed 
as the product of two unipotent matrices, the transition matrix from Jack polynomials to 
monomial polynomials (see [19] Chapter VI). This implies that the span of {m^ : \i\ > n — k} 
again agrees with the span of {d\- : Ai > n — k}, which are rescaled analogues of x X ( see 
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A Sekiguchi-Debiard operator over other bases 

Having seen the probabilistic interpretation of the second order differential operator ffl2|) ex- 
pressed in the power sum symmetric basis p\, it is natural to consider the following question: 
are there other bases of the symmetric polynomials A n over which Lg has natural probabilistic 
interpretation? Here we examine the remaining four fundamental bases: monomial, elemen- 
tary, and complete. The action of L\ on the monomial basis m\ is well-known to be strictly 
upper triangular ([19] p. 317), when the rows and columns of the Markov matrix are arranged 
in a total order compatible with the natural partial order on the set of partitions V n of n: 
// < A if //i + . . . + fi r < Ai + . . . + A r for all r. In particular, if L 2 e does define a Markov 
matrix (meaning the entries are nonnegative), it has a single absorbing state at (1™). 

Next consider the action of L 2 e on e\, the elementary symmetric polynomials. This has 
been studied in detail in [3]. Here we give a quick development avoiding lengthy computa- 
tions. The action of U is easy to describe. For any simple elementary polynomial e r , the 
operators Xj<9j and (xidi) 2 simply collect all the terms in e r that contains the factor Xj. So 
after summing over i G [n], this results in a constant multiple of the identity. Thus to get a 
nontrivial action, one must consider a composite e ri)7 . 2 := e ri e r2 . In this case, one can show 
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that for ri < r 2 , 

n— l 

u ( e n,r 2 ) = 3(1 + n)e rur2 - ^ 2 ( r i + r 2 - 2j)e ri+r2 ^ jd . 

3=0 

Thus C/ is strictly lower triangular with respect to the partial order <. It turns out that the 
action of V on the e\ is diagonal: first of all V satisfies a product rule on e\ = Y[f=l e ^j ! 
by pairing up j ^ k, one also sees that Ve r consists of monomials with no repeated factors, 
hence by symmetry must be a muliple of e r . Thus the following linear combination yields a 
legitimate Markov matrix: 

q 

M e (ci,c 2 ) := c 1 I + c 2 (-U + V). 

Notice that we need to add the muliple of / to make sure that the diagonal entries of M e 
are nonnegative. Also observe that the Jack parameter 9 needs to be non-positive for the 
off-diagonal entries to be nonnegative. It is clear from the description of U and V that this 
Markov chain is absorbing at (n), because the next step either stays in the current state or 
goes to a state corresponding to a partition bigger than that of the current state. 

Next we look at the complete symmetric polynomials h\. First consider the action of 
(X, V) on h r , one of the generators. Since h r = S( r ) is a degree r homogeneous polyno- 
mial, the action of (X,D) = Y^!i=i x A ls si m ply multiplication by r, i.e., any homogeneous 
polynomials are eigenfunctions (X,D). However the operator ((X,D)) 2 acts nontrivially on 
h r : 

n 

^2(xidi) 2 h 4 = -2/ii4 + 10/12,1,1 - 8h 2 2 - 12h 3jl + 28h 4 . 
i=i 

For partitions of more than one parts, the computation gets unwieldy and I haven't tried to 
express U(h ri h r2 ) and V(h ri h r2 ) in terms of h\ explicitly, because of the following numerical 
observation: for A = (3, 2, 1), we have 

Uh\ = 2/i 2 ,i 4 — 8/i 2 2 5 i2 — 2/i 3jl 3 + 14/i 3 2,i + 6/i 3 2 + 6/14,12 + 8/i 4i2 + 10h 5>1 
Vh x = -h 2) i* + 4/i 2 2,i2 + h 3jl 3 + 32/i 3)2 ,i. 

The only linear combination of the above two expressions that yields nonnegative coefficients 
is |[7 + V, which corresponds to 9 = 1. But in that case, the Markov chain is again absorbing 
at (6). So we arrive at the following result: 

Proposition A.l. The operator L 2 Q gives a Markov matrix under the complete symmetric 
polynomial basis h\ for all n if and only if 9 — 1. In that case, the Markov chain never goes 
towards partitions of fewer or equal parts, hence is absorbing at (n). 

Proof. When 9 — 1, h\ are dual to e x with respect to the Jacobi- Trudy identity. Hence the 
walk defined by D\ on h\ can be obtained from the upper-triangular walk on e\ under the 
map A t— > A*; in particular the walk is absorbing at (r). For 9^1, the numerical example 
above suffices to show the associated walk is not positive. □ 

For 9^1, the resulting signed Markov matrix seems to have nontrivial left eigenvec- 
tor corresponding to the eigenvalue 1. I haven't checked if this corresponds to some nice 
stationary distribution on V n ; presumably it will define a signed measure. 
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B Higher order Sekiguchi-Debiard operators 

Throughout this section N will denote the number of underlying variables in the symmetric 
functions and n will denote the weight of partitions, consistent with previous sections. It 
is possible to study higher order differential operators on A^v from the Sekiguchi-Debiard 
operator-valued generating function (see [19J page 328) : 

N N 

D e {X) := asix)- 1 £ e(w)x wS + {w5) j + Bx^) = £ D k 9 X N ~ k . (43) 

weS N j=l k=0 

Below we give complete analysis of D\. It suffices to understand the following two-parameter 
family of operators: 

IN N 

D(\,v;h) = a s (x)- 1 e(w)]J(J2(^)(x j d j r)J2(^d j ) h 

w£S N i=l j=l 3=1 

I N 

3l,— ,31 »=1 i=l 

where A,/i are positive integer compositions and t = £(/i) = £(X). Indeed, it is easy to see 
that (denoting by . . . ,jk) all distinct indices) 



QU " 



D k e =a 5 (xrY: e Em^u^ n 

W£S N (j u ...,j k ) u=0 K ' i=l 



ill 

i=k—u+l 



which can be expressed as a linear combination of D(\, //, /i)'s with |A| + |/i|+/i = k and h ^ 1; 
the factors of the form J2f=i(w5) k or ^2^ = i{xjdj) k all evaluate to constant by symmetry, and 

the operator Y^jLi( x j®j) ac ^ s on P^ by constant multiplication. 

For k = 3, we thus have three operators to consider: D((l), (2); 0), £>((2), (1); 0) and 
.0(0, 0; 3). We compute the action of each on p x below. We need the following three compu- 
tations: 

• a s {x)' l x i d i a 5 {x) = EE^;^%- 

• {xidifpx = Xidi y~] — (xidi)px 3 = y~] — {xidi) 2 px s + V" [{%idi)p>M?idi)P\t\ 

= E x ^ s — + E K\ t x^-^-. 



N <r r+1 1 T r+1 - T r+1 1 r 1 

i=l jjki ■> i^Lj i^j u=0 ri+r 2 =r: 

r,;>l 
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From this we have 

tn\ , v ,._, 



N t{X) , 2 A s -1 

£>((!), (2); 0)p x = a 5 (x)- 1 ^[xidMx^xAfpx = Px(J2 f K 2N ~ A - " X ) + E 



\ \ A s +At — 1 

+ y A f A* [(2iV _ As _ At _ 1) PA £± A i+ v- P,PA 3+ A t - M ]) _ 

It is worth pointing out that the sum Y^=^~ l P "p^p^*~" con tains a constant term (when 
u E {\ s ,\ t }). 

Next we consider D((2), (1); 0). Again we collect some computations below: 
• a 5 (x)" 1 (a;j(9i) 2 a < 5(a;) = (rc^) 2 loga^z) + a; 2 (<9i log a^a;)) 2 

2 2 

_ y^ Xj, y^ 

j^i [ l ]) jjii k^i,j [Xt l Xk) l 3 > 

2 



, :i X i X j j^ k . ( X i x j){ x i x k) 
N r +2 r+2 

\ ^ \ ^ %i y^ y^ \ ^ %i 

i=l jyfc: V X » X ^^ 1 Xk) Tc[N]:ieT j^k: ^ l X ^ )[Xl Xk > 

j,k+i |t|=3 j,fceT\{i} 

=2 s r (T) = 2 E m ^ T ) = 2 E E m ^ T ) 

Tc[iV]: Tc[AT]:Ahr Tc[JV]: Ahr: 

|T|=3 |T|=3 |T|=3 <W<3 

=M N o 3 ) E m * + C\~ 2 ) E "»a+( jv 2 " 1 ) E 

X 7 Ahr: V 7 Ahr: V 7 Ahr: 

£(A)=3 £(A)=2 ^(A)=l 

where s\(T) denotes the Schur polynomial over the variables indexed by T, and simi- 
larly for m\. 



E mA= E E- r ' : -'t /(A ' ^ + 2^N A2/ (Ai = a 

Ahr: Ahr: ijtj 

(A)=2 £(A)=2 

| E (^i^ - Pr) = ^ E 



2 2 

ri+r2=r: ri+r2=r: 
r;>l ri>l 



r _ 1 

Pt\T>T2 ^ Pr 



E mA = ^ E PnPrJPrz - 3 (n - l)p ri p r2 + (r - l)(r - 2)p r ] 

Ahr: ri+r 2 +r 3 =r: n+r 2 =r: 

(A)=3 n>l r t >l 

1 V- r-2 ^ (r-l)(r-2) 
t: 2^ PnPr 2 Pr 3 1^ PriPr2 + fi ' j 



ri+r2+r-3=r: n+r 2 =r: 
ri >l r t >l 
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Putting everything together we have 



2 / / j jriuri-2 ■ g / Vr\PriVrz 

s=l n+r 2 =A s : r 1 +r 2 +r 3 =\ a : 

+ ((AT - l)(iV - A s ) + (2A,-3)(A S -1) K] _ 
Finally we compute the action of -D(0, 0; 3): 

AT £(A) AT £(A) 

d(0,0;3) Pa = Y^w£^— = D^EW^— + E A *W S+At — ]] 

At ^(A) 

tttt Px a % Px s Px t X t Px sP x tP x u 

U Px ° PXt ^ PA ^ Px ° PXtPK 

where Y2(stu) denotes summation over all distinct triples. 

Next we compute the operators D((l), (1); 0) and D(0, 0; 2). 



A" 

D((l),(l);0) Px = a s (x)- 1 ]T e{w)x w5 Y J (w8) i {x i d i )p x 

w£Sm 1=1 

*(A) 



EAs( 1 E ^ + ^, (44) 



s=l ri+r 2 =A s rAs 

rj>l 

JV 

£>(0,0;2)p A = o,(x)- 1 E e(w)x w5 J2(xidi) 2 Px 

weSN »=i 

£(A) 

=p«+E^&. («) 

s=l ^A s fA t 
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We can now compute the action of Dj on power sum polynomials (see [19] Example VI. 3. 3(e)). 

ui i i i i 

i it 

+ g[ lh " U((l), (1); 0)] + y [n 2 - D(0, 0; 2)]p A 

=(^V + ^(iv - 1) + ^-D^-^O^-D j _ eD{{ll 0) _ ^ (0j 0; 2)K 

(46) 

Similarly we can compute D| using the inclusion-exclusion principle: 
DqPx =a s (xy 1 ^2e(w)x wS ^ [-(w5) i (w5) j (w5) k + -(w5) i (w5) j x k d k 

w {i,j,k) 
q2 q3 
+ — (w5)i(Xjdj)(x k d k ) + —{Xidi){Xjdj){x k d k )) 

=(cjv(3, 6) + ^[2,0(0, 0; 3) - 3n£>(0, 0; 2)] + ?-[2D{(l), (2); 0) - 2n£>((l), (1); 0) 
o 3 

- Q D(0, 0; 2)] + 2 1 [D((2), (1); 0) - Q D((l), (1); 0)]K 

where Cjv (3, 0) = J [ (?) 3 - | (*) ) + 2 (?) 2 ] + f [ © \- \ ft) n] + £ (?) n= + ^ . Unfortu- 
nately I cannot extract any natural interpretation of Markov chains from the right hand side. 
This is not so surprising since the composite walk Pq for k > 2 does not correspond to some 
affine transformation of the Metropolization of Pf with respect to the measure MED(0). 
Nonetheless this gives a new Markov chain that converges to the Multivariate Ewens distri- 
bution with parameter 6' -1 , since the operators D r e are simultaneously diagonalized and the 
left eigenfunction corresponding to the eigenvalue 1 is simply the stationary distribution. 

I also computed a numerical example using the symmetric reduction function in mathe- 
matica and the SF package in maple. We take the power sum polynomial p\ with A = (3, l 2 ). 

Dfazpl = -39p 2 p\ + (339 + 35 + 79 2 )p 3 p\ - 69 2 Pm - p 2 p 3 9 2 

D 3 e p 3 p 2 = (2/3)9p\ + (-A9 2 - W)p 2 p\ + (229 2 + (73/6)# 3 + 50 + (307/3)% 3 P? + ^PipI 
+ (49 3 - 209 2 )p m + (-(4/3)# 3 - 29 2 )p 3 p 2 + 69 3 p 5 

Dg o D 2 B p 3 p\ = 39 2 p\ + (-219 - 49 s - \§9 2 )p 2 p\ + (505# 3 + 1225 + 23100 + A9 A + lh7% 2 )p 3 p\ 
+ 2A9 3 p 1 p\ + (-410 3 - Q9 4 - A29 2 )p 4Pl + (-Q9 3 - 9 A - 19 2 )p 3 p 2 + 30/4 

This example shows that D 2 ,, Dg, Dj, o D 2 e and id are independent operators on A^v- Notice 
also that DqPx has positivity issues: the partitions of Cayley distance 2 from the starting 
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partition A are always positive, whereas the ones that differ from A by one transposition 
might become negative. So in order to make Dg into a Markov matrix, one needs to add 
a sufficiently negative multiple of Dg. We have not tried to compute the optimal multiple 
here since we are unable to glean any nice pattern from the numerical example above; in 
particular, the coefficients cannot be made into simple powers of 9. Observe that Dg o D 2 e p\ 
also has positivity problem, but it is much easier to fix, since one can simply add a multiple 
of the identity to Dg as in the case treated by [5]. 



C Compositions of D 2 e for different 6 values 

In general the eigenvalues and eigenfunctions of a Markov chain can be highly intractible, 
due to the need to solve for high degree polynomials. For instance, the Metropolis chain 
based on 3-cycle shuffle on V n already requires taking square roots for n — 4. 



Mf\9) 



( 



e 2 
V o 



o 

1/8 

1 


9 2 




3/4 






1 




1-9 2 





1/8 



1-9 2 



\ 



The eigenvalues are 



1,1, 



-6 2 , — ( 1 - 86 2 - V193 - 208fl 2 + 649 A ) , — ( 1 - 86 2 + Vl93 - 2089 2 + 640^ 
16 V / 16 V j 

The following result was discovered in numerical experiments: 

Proposition C.l. For any Laurent polynomial p in m variables, p(Dg , . . . , Dg ) gives rise 
to a Markov chain on the set of partitions V n , with eigenvalues, and left and right eigenvectors 
given by rational functions of 9\, . . . , 9 m . In particular, the stationary distribution is given 
by rational functions of 9i 's also. 

Proof. When expressed in the monomial symmetric basis, D 2 , is unipotent (upper triangular 
with l's on the diagonal), with respect to any total ordering on V n compatible with the 
natural partial ordering <, whereby A < n if Ai + . . . + A r < \i\ + . . . + /v for all r (see [TS] 
page 317 equation (3.7)). Thus fixing this total-ordering, any Laurent polynomial of Dg.'s is 
clearly still unipotent. The eigenvalues are simply the diagonal entries, and the eigenvectors 
can be computed using simple row reduction, which also results in rational components. □ 

The above result is clearly also true for Dg in general and even Macdonald operators. 
Thus in principle, one can compute the C 2 mixing time of Markov chain of the form Dg ± Dg 2 , 
whose stationary distribution can be quite complicated: for n = 3, the stationary probabili- 
ties are 



{ 



9 2 (0 X (3 - A9 2 ) + 9\{-l + 9 2 ) + 9 2 ) 



3 (-(-3 + 9 2 )9 2 + 9 X (1 - 4# 2 + 9 2 )) 



1 + 9fl 2 + 9 2 ) 



8 + 9 2 + 9 2 (-l + 9 2 )9 2 - 9 1 

-2fl 1 + 2(-2 + fl 2 ) 2 

8 + 9 2 + 9 2 (-l + 9 2 )9 2 - 9x (-1 + 9# 2 + 9 2 ) 



8 + 9 2 + 9l(-l + 9 2 )9 2 - 9 1 



■1 + 9fl 2 + 9 2 ) 



} 
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D Extensions to other root systems 



The appropriate generalization of the Laplace-Beltrami operator to root systems other than 
A N is given by the Heckman-Opdam operator (see [13] and [3]): 

N 

L n (k, R) = A + kV n := ^2 9 l + E Ka cot H a / 2 ) d a- 

i=l a£i?+ 

where R denotes an arbitrary root system, R + a designated set of positive roots, and K a is 
called a multiplication function, invariant under the action of the Weyl group on R + . For 
more on root systems and Weyl groups, consult the first 3 chapters of [UJ as well as chapter 
19 and 20 of @]. 

We study root system Dn in detail here. These come from the irreducible decomposition 
of the adjoint representation of the maximal torus in the compact Lie groups SO(2N, R). 
The positive roots can be chosen as the set {x iXj j Ob j • 1 < ^ < j < N} on the maximal 
torus; in the associated Cartan subalgebra (the Lie subalgebra corresponding to the maximal 
torus), they become {U —tj,ti + tj : 1 < i < j < N}. The appropriate analogue of the power 
sum polynomials appears to be the power sum symmetric Laurent polynomials given below: 

N N 
i=l i=l 

where Xi = e*\ And as in the case of A N , p x = Yli=i \p\v By direct computation we have 

a-1 

coth(0/ '2) d a p a = 2a ^ E cosh (^) cosh((a - £)tj) 

0&R+ ijtj £=0 

a— 1 a— 1 

= (2an - a 2 - a)p a + 2a 2j PtPa-l ~ a E P a ~ 2 ^ 

i=i e=i 



a-i 



EJ 2 Xj=i 1)/2 Pa-2£ if a is odd 

Pa-2£ — S v^(a-2)/2 . f . 

N + 2 2_Je=i Pa-2t it a is even 



£=i 

LfJ 

= 2^Pa-2i, 

if we define po := iV/2. 

Therefore for n = |A| = A«, 

*(A) Ai-l LA»/2J 

E coth(a/2)9 aPA = ((2iV - l)n - ^ \ 2 )?a + £ — 2A ? ,[^ PO»*-t ~ E 
aei?.+ i=l i=l £=1 €=1 
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Restricting to partitions of the top grading, n, clearly the transition coefficients are affine 
transformation of those in the case, hence nothing new is obtained this way. There are 
several pathological features regarding the action of L^(D, k) := A + Yl a &R + K <* coth(a/2)9 Q; 
on the power sum analogues of symmetric Laurent polynomials (see the toy example below): 

1. There are no easy ways to make the entries all positive, 

2. The row sums are not the same. 

Thus it remains difficult to interpret the full transition matrix as a Markov kernel. For root 
system D N , there is only one Weyl orbit, hence K a = k. The Heckman-Opdam functions 
have rational transition coefficients to this power sum Laurent basis, as illustrated by the 
following numerical example (Vk denotes the set of partitions of k): 



M K (N)\ v . iUVl := 

/ 9 + k(-12 + 6iV) 12k -6k \ 

2 5 + /<-8 + 6iV) Ak -2-2kN 

3 3 + K(-6 + 6iV) -3 

\ 1 + k(-2 + 2N) J 

where the columns and rows are indexed by (3), (21), (l 3 ), (1). The eigenvalues are very 
clean: 



3 + 6k(-2 + N), 9 + 6k(-1 + N), 5 + «(-8 + 6iV), 1 + 2«(-l + N) 

The left eigenvectors are rational functions of k, which we display as rows of the following 
matrix: 



/ 



-5+6+2N 
3(-l+JV) 



(26>(-l+26>+7V) 

3(l+26»+V) 
6>(-3+26>+2V) 
0+26 2 -2N+20N 



-5+6+2N 
-1+N 



-2(-l+20+AT) 
1+29+N 
(2(-l+g)(-3+2g+2V) 
6+2d 2 -2N+20N) 



/ 2(-5+0+2N) 
V 3(-l+JV) 


(4(-l+26>+V) 
36»(l+26»+V) 
12-8e~8V 
6»+2e 2 -2V+26»V 



Here 



corresponds to the parameter in the case. 



We have also tried to adjust the diagonal entries to make the row sum equal to 1; the 
resulting matrix however does not have rational eigenvalues in the entries. 
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